Analog gravity from Bose— Einstein condensates 
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We analyze prospects for the use of Bose-Einstein condensates as condensed-matter systems 
suitable for generating a generic "effective metric", and fo r mimicking kinematic aspects of general 



relativity. We extend the analysis due to Garay et al, |gr-qc/000201^ , ^r-qc/000513ll] . Taking 
a long term view, we ask what the ultimate limits of such a system might be. To this end, we 
consider a very general version of the nonlinear Schrodinger equation (with a 3-tensor position- 
dependent mass and arbitrary nonlinearity) . Such equations can be used for example in discussing 
Bose-Einstein condensates in heterogeneous and highly nonlinear systems. We demonstrate that 
at low momenta linearized excitations of the phase of the condensate wavefunction obey a (3-1-1)- 
dimensional d'Alembertian equation coupling to a (3-|-l)-dimensional Lorentzian-signature "effective 
metric" that is generic, and depends algebraically on the background field. Thus at low momenta this 
system serves as an analog for the curved spacetime of general relativity. In contrast at high momenta 
we demonstrate how to use the eikonal approximation to extract a well-controlled Bogoliubov- 
like dispersion relation, and (perhaps unexpectedly) recover non-relativistic Newtonian physics at 
high momenta. Bose-Einstein condensates appear to be an extremely promising analog system for 
probing kinematic aspects of general relativity. 
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I. INTRODUCTION 

Progress in understanding classical general relativity (GR) and quantum field theory in curved spacetime (not to 
mention quantum gravity itself), suffers greatly from the lack of experimental feedback. Moreover, direct experimental 
probes of many important aspects of these theories (such as, for instance. Hawking radiation from a black hole) appear 
to be very far from current (and even foreseeable) technologies. In this regard, the possibility of using condensed 
matter systems (such as, for example, Bose-Einstein condensates [|l|, |^) to mimic certain aspects of GR could prove 
to be very important. The basic idea of a condensed matter analog model of GR is that the modifications to the 
propagation of a field/ wave due to curved spacetime can be reproduced (at least partially) by an analog field/wave 
propagating in some material background with space and time dependent properties. 

The concept of a condensed matter analog model was first explored in a little-known paper by Gordon, where he 
worked within the context of the optical properties of dielectrics After lying fallow for a considerable period, 
during the last decade or so this idea has been revived and elaborated, mainly in the context of considering the 
propagation of sound waves in a movin gfluid 11 I I, 0, |, |, |ll|, H. Several other analog systems have 
also been analyzed; see for example l]!^, p^, and also the mini-review by Jacobson The search for new analog 
models, and exploration of known analog models, is ongoing jl^, Among these models, a particularly 

promising one that will be the focus of this article is the recent proposal of Garay, Anglin, Cirac, and ZoUer based on 
Bose-Einstein condensates (BEG) ^. 

Over the last few years a remarkable series of experiments on vapors of rubidium |2l[, lithium |^^ , and sodium 
psf have led to a renewed interest in the phenomenon of Bose-Einstein condensation ]24| , p5| . In these experiments 
gases of alkali atoms were confined in magnetic traps and cooled down to extremely low temperatures (of the order 
of fractions of microkelvins) . In order to observe the BEG the whole gas was allowed to expand, by switching off the 
trapping potential, and monitored via time-of-flight measurements made with optical methods. The signature of the 
BEG was a sharp peak in the velocity distribution for temperatures below some critical value (see for an extensive 
review on the subject). 

As Garay et al have shown, perturbations in the phase of the condensate wavefunction satisfy, in the low-momentum 
regime, an equation equivalent to that of a massless scalar field in a curved spacetime (the d'Alembertian equation 
A(f) = 0), but with the spacetime metric replaced by an effective metric that depends on the characteristics of the 
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background condensate. Present-day experimental achievements, and the rapid development in magnetic trapping 
techniques, seem to illuminate a viable path to experimentally reproducing important general relativistic features 
such as ergoregions and event horizons , D . 

In this paper we wish to explore the Bose-Einstein system in more detail, formally extending the Garay et al analysis 
as much as possible. We will analyze a number of physically conceivable extensions of the usual Gross-Pitaevskii 
approximation. In particular, we will consider arbitrary non-linear interactions and anisotropic mass-tensors, both 
depending explicitly on (space and time) position. In this scenario we will investigate how close we can come to 
mimicking a "generic" gravitational field. (The simple analog models based on ordinary fiuid dynamics are somewhat 
limited in this regard because the spatial slices are always conformally flat ^, |l^). After including these 
generalizations the only significant constraint left on the effective metric is due to the irrotational nature of the 
condensate 3-velocity. Within this scenario we will then show how various quantum corrections come into play in the 
whole analysis, distinguishing three useful regimes: the quasi-classical, low-momentum, and high-momentum regimes. 

By using eikonal techniques we shall investigate the high-momentum dispersion relation of the collective excitations 
of the condensate. In particular we shall study it with regard to the possibility of analyzing the "Hawking radiation" 
that may be emitted from "horizons" that form in this system. We find a Bogoliubov-like dispersion relation of 
the superluminal (more properly, supersonic) type. For situations with both an outer and an inner horizon, this 
dispersion relation seems to led to exponential amplification of the radiation flux coming from the outer horizon , 
(a so-called black hole laser), in agreement with the unstable behaviour of some of the solutions found by Garay et 
al. (They also found some metastable solutions whose origin might be related to the speciflc boundary conditions of 
the configurations considered.) Finally, we discuss the way in which in the high-momentum limit one recovers the 
underlying non-relativistic Newtonian physics. 

II. BOSE-EINSTEIN CONDENSATES 

In a quantum system of A'' interacting bosons the crucial feature of a Bose-Einstein condensate (BEG) is that it 
corresponds to a configuration in which most of the bosons lie in the same single-particle quantum state. This system 
can be suitably described in a second quantization framework by a many-body Hamiltonian of the form |^ : 

/r ^- 2 1 1 /" 

dx^Ut,x) V^ + Vo^tM $(t,x) + - /dxdx' $^(t,x) $^(t,x') l^(x-x') $(i,x') *(i,x). (1) 
2m J 2 J 

Here T4xt(x) is some confining external potential, y(x — x ) is the interatomic two-body potential (other possible 
multibody interactions are neglected at this stage), and m is the mass of the bosons undergoing condensation (in 
current experiments these bosons are actually alkali atoms). Finally, ^'(t,x) is the boson field operator. 

Although the quantum state for the Bose-Einstein condensate, as well as its thermodynamic properties, can in 
principle be exactly computed from JH ) , it is clear that for large ensemble of atoms this approach can become 
impractical. It was Bogoliubov |28| p9| who first recognized that a natural ansatz for studying such a system is 
the mean-field approach, which consists of separating the bosonic field operator '${t,x.) into a classical condensate 
contribution ?/'(t,x) plus excitations ^(t,x): 

$(t,x) =7A(i,x) + (?(t,x). (2) 

Here ip{t,x) is defined as the expectation value of the field ip = (^(t,x)). It is sometimes referred in the literature 
as the "wave function of the Bose-Einstein condensate" . Its modulus fixes the particle density of the condensate, 
p — N/V, in such a way that \tp{t,x.)\'^ = p(t,x). 

The Bogoliubov decomposition (H) is particularly useful when the number of atoms which do not lie in the ground 
state of the condensate is small. In this case it is in fact possible to consider the "zeroth-order approximation" to 
the Heisenberg equation given by the many-body Hamiltonian (|l|) by just replacing ^(t, x) with the condensate wave 
function 



d \ 



2m 



+ Kxt(x) + J dx' V^(x', i) F(x' - x) V(i, x') 



V'(i,x). (3) 



The next approximation is to assume the intermolecular interaction term V{x' — x) represents a short range inter- 
actions. (In helium-based Bose-Einstein systems this is not the case and considerably more complicated analysis 
is required.) Part of the theoretical interest in the previously cited heavy-alkali-atom Bose-Einstein condensates is 
that these are extremely dilute condensates; systems in which at low energy only binary collisions are relevant. This 



3 



permits one to model the interaction with a short-distance delta-Uke term times a unique self-couphng constant A 
which is determined by the s-wave scattering length a and the atomic mass m 



V^(x'-x) « A(5(x'-x), (4) 
A = (5) 



The use of the approximate potential in equation leads to a standard closed-form equation for the weakly 
interacting boson condensate 

i?i^V'(i,x) - (^-|^V2 + Fext(x) + A |^(t,x)|'^ VXi,x). (6) 

This equation, most commonly known as Gross-Pitaevskii (GP) equation, (sometimes called the nonlinear Schrodinger 
equation, or even the time-dependent Landau-Ginzburg equation), can be associated to an effective action of the form 

y"did3x|^*(^i;iat-f |-v2-Kxt(x)^V^-iA |V'(i,x)|^|, (7) 

which takes the name of the time-dependent Landau-Ginzburg action. (There is some disagreement on terminology 
here: some authors prefer to use the phrase "time dependent Landau-Ginzburg equation" for equation (^) only after 
formally discarding the i. This has the net effect of replacing the Schrodinger differential operator by a diffusion 
operator. We always keep the i; all our physics is oscillatory rather than diffusion-based.) 

The Gross-Pitaevskii equation describes, in a simple and compact form, the relevant phenomena associated with 
BEG. In particular, it reproduces typical properties exhibited by superfluid systems, like the propagation of collective 
excitations and the interference effects originating from the phase of the condensate wave function. It is precisely 
from this (non-relativistic) theory that we shall now see how an analogue of GR can be constructed, and possibly 
used for experimental purposes. 



III. GENERALIZED GROSS-PITAEVSKII EQUATION 



The aim of our investigation is to explore the ultimate extent of the ability of BEG systems to mimic general 
relativistic ones. In order to do this in the most general way we shall now introduce a generalization of the nonlinear 
Schrodinger equation (|^). In particular, we shall formally consider a series of generalizations each of which is in 
principle allowed for systems undergoing Bose-Einstein condensation, though they are not (yet) commonly encountered 
in the experimental literature. 

— (1) The first generalization we will make is to replace the quartic iA|i/)(t, x)!"* by an arbitrary nonlinearity 
Tr{ip*ilj) = tt{\iP\'^). We note in particular that for two-dimensional systems Kolomeisky et al have argued that 
in many experimentally interesting cases the nonlinearity will be cubic or even logarithmic in ji/'P- (Although in 
two-dimensional systems standard Bose-Einstein condensation does not occur it is experimentally almost certain that 
"quasi-condensates" exist, and possess collective excitations treatable in an analogous way |3l|.) 

— (2) The second generalization is that we will also permit the nonlinearity function to be explicitly space and 
time dependent: tt — > 7r(x, [V'*V']) with x = (t,x). While intrinsic properties of the bosons that we are trying to 
condense (like the scattering length and mass) will typically not change from point to point, we can certainly envisage 
a more general situation in which the condensate is either selectively "doped" , or perhaps (as suggested in Garay et 
al 0, ^) physically constrained to move thorough a narrow tube of varying cross-section, thereby altering its effective 
properties in a position-dependent way. (We do not mean to imply that this would be easy, or that such techniques 
are "just around the corner"; instead we are interested in seeing how far we might ultimately be able to push this 
system.) 

— (3) The third generalization we will make is to permit the mass to be a 3-tensor: m — s- rriij. Such anisotropic 
masses are well-known from condensed matter physics where they are most typically encountered in effective mass 
calculations for electrons immersed in a band structure (see, for example p^). They have also been discussed for the 
case of excitons (electron-hole couples held together by Goulomb attraction) in BEG for semiconductors. The doping 
structure of the semiconductor and its anisotropics would give place to an effective mass matrix for the paraexcitons 
(singlet excitons) at least in the low momentum approximation Q . Formally we shall consider the possibility for 
anisotropic masses in the more general context of the nonlinear Schrodinger equation, regardless of whether or not we 
are dealing with a band structure or even a Bose-Einstein condensate. (If you wish, you might like to think of this 
as a liquid-crystal BEG.) 
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— (4) The fourth generahzation we will make is to also permit the 3-tensor mass to depend on position (both time 
and space). Again we might like to think of a doping gradient or similar situation. Mathematically this has the effect 
that the mass matrix must be viewed as a "metric" on a curved 3-dimensional space. It is convenient to introduce an 
arbitrary but fixed (time and space independent) scale fi with the dimensions of mass and then write 

= ^I'-^'^h.j (8) 

with ^^"^hij being a properly dimensionless 3-metric. Note that the introduction of /i is a mathematical convenience, 
not a physical necessity, and that all properly formulated physical questions will be independent of fi. 

— (5) The fifth and last generalization will be to allow time-dependence for the confining potential 14xt = Vcxtit, x). 
(This is already implicit in the analysis of Garay et al, but is somewhat non-standard from the usual condensed-matter 
viewpoint.) 



We do not claim that any of these generalizations will be easy to achieve experimentally, it is sufficient for our 
purposes that they are at least physically conceivable. After this series of generalizations we obtain an action which 
now reads: 



S 



2fi 



(3) 



Rih)-V,^t{t,^) U-^(a;,|Vn 



Here is the 3-dimensional Laplacian defined by 



det[(3)/i] 



/det[(3)/i] 

where [(3)/j-i]*i jg the inverse of the 3-metric '■^^hij. Additionally note the presence of the DeWitt term 



(9) 



(10) 



(3) 

2/i 



R{h) 



(11) 



involving the dimensionless parameter ^ and the 3-dimensional Ricci scalar — this term arises from operator-ordering 
ambiguities in going from the "flat space" metric to "curved space" (going from position-independent m to a position- 
varying effective mass) |35, 3^. We include the DeWitt term here for completeness, and because it should be included 
as a matter of principle, but note that it is unlikely to lead to experimentally measurable effects. 

Varying this generalized action with respect to ip* now gives the "generalized" nonlinear Schrodinger equation that 
will be of central interest in this article: 



Ahi'it, x) - 1^ '^^^R{h) i,{t, x) + Fext(i, x) ^(t, x) + ^'{r^) m x) 



(12) 



We will now demonstrate that this non-relativistic generalized nonlinear Schrodinger equation has a (3-fl)-dimensional 
"effective" Lorentzian spacetime metric hiding inside it. 



IV. MADELUNG REPRESENTATION AND THE HYDRODYNAMIC LIMIT 

The so-called Madelung representation ^ of a Schrodinger wavefunction consists of writing 

iPit, x) = ^/p{t,^) exp[-i6i(t, x)/^]. (13) 

The factor of h is introduced for future convenience (it suppresses h as much as possible in the following equations). 
This implies that 6 has the dimensions of an action. The Madelung representation is well-known in the context of the 
ordinary linear Schrodinger equation, and generalizes to the present situation without difficulty. 

Garay et al |l|, ^ substitute the Madelung representation into the Gross-Pitaevskii equation, and take the real and 
imaginary parts. We could do the same thing here with the generalized nonlinear Schrodinger equation. Alternatively 
we could insert the Madelung representation directly into the action, and vary with respect to and p to deduce 
Euler-Lagrange equations. Either way, you will get two equations: 
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(1) Continuity: 



dtp+-V-{pV9)^0. (14) 



Here and hereafter the V denotes a covariant derivative with respect to the 3-nietric /i^ . If we define a "velocity" 
(and at this stage this is a purely formal step) 

(v)' = = [m-'Y^ V,0, (15) 

M 



then this "velocity" is actually independent of /i, and equation ( [14D above is formally equivalent to the usual equation 
of continuity in a curved 3-space 

^p + V-(pv)=0. (16) 
— (2) Quantum analogue of the Hamilton-Jacobi equation: 

l.(Vf?)2 + y,^, + ,'_5! f^+e (3) !?(/,)')= 0. (17) 



dt 2fj, 2fi \ y/p 

Here and hereafter (V0)^ denotes the 3-metric scalar inner product (3)[/j-i]y Q.g Q^Q^ 
The quantity 

is a generalization (because it now includes the 3-dimensional Ricci scalar term) of what is often called the "quantum 



potential" 1 42 , 4^ . (Note that the quantum potential is actually independent of /i because rescaling p simultaneously 
rescales hij in a compensating manner.) With our conventions this is now the only place where h appears. 
In terms of the the velocity field ( p^ ) we can write the Hamilton-Jacobi equation in the form 

I K-(v)^] + V, (^HhnS^p^ + v,^, + n' + Fq) = 0. (19) 

(Note that this form makes the /^-independence of the physics manifest.) The two real equations (|l^) and (|r^ [or 
alternatively equations (^6|) and ( | l9 [) i subject to the definition (^5|)] are completely equivalent to the generalized 
nonlinear Schrodinger equation ([1^)7 

An interesting physical regime is that in which one can safely neglect the quantum potential term. This approxi- 
mation can be justified either as the classical limit of the theory (it corresponds to neglecting all terms with powers 
of h) or as the regime of strong repulsive interaction among atoms. In the latter case the density profile can be safely 
considered smooth and hence it is reasonable to neglect the kinetic-pressure term A^y^/y/p. This, together with the 
smallness of the DeWitt term, permits one to discard Vq . 

In any case one can see that neglecting Vq, the equations ( p^ ) and ( p9| ) have the form typical of those for supcrfluids 
in the T — > limit. In particular we can see, by the absence of any term proportional to v x (V x v) in equation ([is|), 
that the equations we are working with are automatically vorticity free. This is generally a necessary assumption 
in order to obtain tractable equations in the analog of GR from standard hydrodynamics |ll[; here it is a free 
byproduct of the GP equation (and this conclusion is does not depend on the assumption of neglecting Vq). 

The hydrodynamical form of the equations ( |l6|) and (|l^) allows to describe the Bose-Einstein condensate as a gas 
whose pressure and density are related by the barotropic equation of state 

Pip)^7T'ip) p-7r{p). (20) 

Therefore it is also possible to formally define a local speed of sound by varying this pressure with respect to the mass 
density of the condensate "fiuid" g ^ pp 

(21) 

dg n 
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We say formally because for the general anisotropic case this "velocity of sound" is not physical. Physically, there 
will be three principal sound velocities in three orthogonal directions, and, as we have mentioned before and will see 
again, the formally convenient parameter does not appear in any true physical result. This does not diminish the 
convenience of introducing the parameter fi for intermediate stages of the calculation. 

It is also interesting to give a more quantitative estimate of the magnitude of this sound velocity for the Bose-Einstein 
systems based on alkali atoms. In particular we can consider the specific case of the rubidium gas with trivial mass 
tensor = m) and standard interaction term 7r(t, x, = 7i'(t, x, p) = ^Ap^, so that 7r'(t, x, = A|V'(i, x)p = Xp. 



[Here A is given by equation (||).] In this case the equation of state becomes P 
the well-known form 



and the speed of sound takes 



2h 



(22) 



For the rubidium gas one has a( Rb) « 5.77 nm, ■m{° Kb) « 86.9 u and in standard BEC experiments p ~ 
10^^ cm""^ These numbers lead to a value of the speed of sound c ~ 6.2 x 10^"^ m/s « 6 mm/s. This is indeed 
one of the lowest speeds of sound one can experimentally obtain. We shall see in what follows how this number can 
play an important role in the simulation of gravitational phenomena in BEC systems. 



FLUCTUATIONS 



Now that we have seen how it is possible, at least in some appropriate regime, to introduce a hydrodynamical 
interpretation of the condensate equations, and how a speed of sound can be meaningfully introduced, we are naturally 
lead to investigate the propagation of fluctuations in the condensate. 

In order to purse such an investigation we shall linearize the equations of motion ( p^ ) and ( p^ ) around some assumed 
background {po,6q). In particular we shall set p = po + epi + O(e^) and 9 — 6o + e^i + 0{e^). Then, we will be 
left with two equations for the background configuration plus two more (often called Bogoliubov equations) for the 
linearized quantities. Linearizing the continuity equation results in the pair of equations 



dtpo + -V • (po V0o) = 0, 
p 

dtpi + -V • (pi VOq + Po VOi) = 0. 
M 

Here and hereafter all inner products (a • b) are all calculated using the 3-metric (/ly a* V). 
Linearizing the Hamilton-Jacobi equation we obtain the pair 

dtOo + ^(V^o)' + Kxt + + ^q(po) = 0, 



(23) 
(24) 

(25) 



1 T> 
dt9i + -VOo ■ Ve, + 7r"(po) — ^2Pi = 0. 
p Ip, 



(26) 



Here D2 represents a relatively messy second-order differential operator obtained from linearizing the quantum po- 
tential, explicitly: 



n — -'- -3/2 FA / +1/2m I 1 -1/2 A / -1/2 n 

D2P1 = --Po [^^h[Po )\Pi + ^Po ^hiPa Pi)- 
The linearized Hamilton-Jacobi equation may be rearranged to yield 



(27) 



Pi 



^"{Po) - ^ ^2 



(28) 



Now substitute this consequence of the linearized Hamilton-Jacobi equation back into the linearized equation of 
continuity. We obtain, up to an overall sign, the wave-like equation: 



P' 



^"{Po) - ^ ^2 



dtei 



1, 
p 



^"{Po) - ^ ^2 




V0n ■ ve 




(29) 
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This wave-like equation describes the propagation of the hnearized Schrodinger phase 0i. [The coefficients of this 
wave-hke equation depend on the background field {po,Oq) that you are linearizing around.] Once 9i is determined, 
then equation determines pi. Thus this wave equation completely determines the propagation of linearized 
disturbances. The background fields po and Oq, which appear as time-dependent and position-dependent coefficients 
in this wave equation, are constrained to solve our generalized nonlinear Schrodinger equation. Apart from this 
constraint, they are otherwise permitted to have arbitrary temporal and spatial dependencies. To simplify things 
construct the symmetric 4x4 matrix 



j-iO 



(30) 



(Greek indices run from 0-3, while Roman indices run from 1-3.) Then, introducing (3 -I- l)-dimensional space-time 



coordinates 



Here 



{t]x^) — the above wave equation ( |29| ) is easily rewritten as 

df.if^" a,0i) = 0. 



r 



7r"(po) 



— D2 
2p 



2p 



A; Co 



p 

-2 -. -1 



V'(,o) - 2^ D2 



Po 



p- 



p 



P 



(31) 

(32) 
(33) 
(34) 
(35) 



Thus Z^" is a 4 X 4 matrix of differential operators (the differential operators in question consistently operating on 
everything to the right). Note that the precise placement of the h^-' above is immaterial since the operator D2 is 
built using V, the covariant derivative associated with h^^ . This remarkably compact formulation ( |3l| ) is completely 
equivalent to equation ( |2^ ) and is a much more promising stepping-stone for further manipulations. 

The major obstruction to interpreting this wave equation in terms of Lorentzian geometry is the fact that /^"^ is 
itself still a matrix of differential operators, not functions. We now consider several different approximations (valid 
in different regimes) which have the effect of replacing these differential operators by functions. After making those 
approximations, the remaining steps are a straightforward application of the techniques of curved space (3 -I- 1)- 
dimensional Lorentzian geometry. (See, for instance & uM ) 



VI. QUASI-CLASSICAL APPROXIMATION 



The most straightforward approximation we could make is to simply neglect D2 completely: This is actually what 
is done in the analysis of Garay et al |^. We (and they) justify this approximation by pointing out that D2 is 
always multiplied by ti^ and in this sense is suitably "small" . We shall call this the quasi-classical approximation. 
The wave equation then simplifies to 

f {dtei + [i/p]V0o-vOi) \^ ^ 1 ^ f xpo ve, - vop (dtO, + [i/p]V9o ■ ve,) \^ ^ ^ ^^^^ 



TT J P \ 

To simplify things algebraically, we can define, in analogy with equation (|2^), a speed c as 

2 _ Pa 



P 



and the background "velocity" Vq by 



(vo) 



P 



(37) 



(38) 
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(Remember that c is not necessarily the physical speed of sound; it is simply a convenient parameterization; in contrast 
vg really is /i- independent and physical.) Now construct the symmetric 4x4 matrix 



(c2 h'^ -v'avi) 



(39) 



This is now just a 4 x 4 matrix of numbers. 

In any Lorentzian (that is, pseudo-Riemannian) manifold the curved space scalar d'Alembertian is given in terms 
of the (3+l)-metric g^„{t,x) by 



A9 = 



1 



5m {V^ 5^" da 



(40) 



The (3+l)-dimensional inverse metric, g^'^{t, x), is pointwise the matrix inverse of g^v{t, x), while g = det((7^^). Thus 
one can rewrite the physically derived wave equation (pq) in terms of the d'Alembertian provided one identifies 



V^9 9'"' = Fr 

This implies, on the one hand 

det(/'^-) = {^gf = g. 
On the other hand, from the explicit expression (|30|), expanding the determinant in minors 
det(/-^) = [n'T' [(-1) • (c^ - ^^0^) - (-^^o)^] • [c'] ■ [c'] = -cV{n"r 

Thus 

g^^cVin'T; ^g = c^l[j:"f. 
We can therefore pick off the coefficients of the inverse "condensate metric" 



5^''(t,x) 



-1 



(c2 h'^ 



VnV, 



POC 



-1 



(c2 - v'oV^g) 



(41) 
(42) 
(43) 
(44) 

(45) 



This is the effective Lorentzian metric seen by the perturbations of the phase of the condensate wave function. At 
this point let us compare this metric with the acoustic metric of ||, H, |^, |8[ |[ ^ ^ |l2j : 

• (1) The two metrics (acoustic and condensate) possess the same conformal factor, (up to a physically irrelevant 
constant rescaling). In view of this we will just call it the acoustic metric from now on, keeping in the back of 
our minds that the relevant "acoustics" is now the propagation of oscillations in the phase of the condensate 
wavefunction. 

• (2) Note that Garay et al did not keep track of the conformal factor, as it was not needed for the points they 
wanted to make. 

• (3) There are now slightly different physical interpretations for c and vq. 

• (4) There is already a non-flat 3-metric hij present in the analysis, even before the linearization procedure is 
carried out. It is this feature that departs furthest from the previous implementations of the notion of "acoustic 
metric" . 

We could now determine the metric itself simply by inverting this 4x4 matrix. On the other hand, as is by now 
standard, it is even easier to recognize that one has in front of one an example of the Arnowitt-Deser-Misner split of 
a (3 -|- l)-dimensional Lorentzian spacetime metric into space + time, more commonly used in discussing initial value 
data in Einstein's theory of gravity — general relativity. The (direct) acoustic metric is easily read off as 
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Equivalently, the "acoustic interval" can be expressed as 

ds^ ^ dx^" dx" = ^ -c^dt"^ + h,j (dx' - vif dt) {dx^ - vi dt) . (47) 

A few brief comments should be made before proceeding: 

• Observe that the signature of this metric is indeed (—,+,+,+), as it should be to be regarded as Lorentzian. 
There is an interesting physical subtlety here: Some alkali atomic gases have a negative scattering length, that 
is, there are attractive forces between atoms physically leading to the collapse of the BEG. A negative scattering 
length is formally equivalent to an imaginary speed of sound, and in terms of the effective metric is equivalent to 
a Euclidean-signature metric. That is: negative scattering length corresponds to an elliptic differential operator, 
instead of the more usual hyperbolic differential operator. In terms of general relativity, manipulating the sign 
of the scattering length corresponds to building an analog for a signature-changing spacetime. 

• It should be emphasized that there are (at least) three distinct metrics relevant to the current discussion: 

— The physical spacetime metric is just the usual flat metric of Minkowski space 

Vt^v = (diag[-cfight, 1, 1, 1])^^. (48) 

(Here Qight ~ speed of light.) The quantum field couples only to the physical metric ?7^t^. In fact the 
quantum field is completely non-relativistic — \ \vo\ \ ^ Qight- 

— The "mass metric" rriij = fj, hij describing the position-dependent effective mass for the fundamental 
particles described by the generalized nonlinear Schrodinger equation. 

— Fluctuations in the condensate field on the other hand, do not "see" the physical metric at all. Perturbations 
couple only to the acoustic metric g^i, . 

— As is common to all two-metric theories, one could also construct a distinct "associated metric" by the 
prescription 

(49) 

There seems to not be any clean physical interpretation for this object. 

• The geometry determined by the acoustic metric inherits some key properties (such as, for example, stable 
causality) from the existence of the underlying flat physical metric. (See [^j, pl|.) 

• This acoustic metric is now sufficiently general to be able to mimic a wide class of "generic" gravitational fields — 
the presence of the position-dependent 3-tensor mass matrix is crucial to this observation. In Garay et al |^ the 
mass was both isotropic and position- independent, consequently the spatial slices of their analog spacetimes were 
always conformally flat (the same phenomenon occurs in the acoustic geometries of § |, |, |, ||, |l|, H ^ ) . 
This is the fundamental reason we have gone to the technical trouble of adding a position-dependent 3-tensor 
mass. The only significant restriction on our version of the effective metric is that the three- velocity is irrotational 
(zero vorticity; curl-free). 

• The major differences with the Garay et al analysis is that their nonlinearity was strictly quartic, their mass 
both isotropic and position-independent, and that they further approximated the current quasi-classical approx- 
imation by going to a "geometrical optics" version thereof that permitted them to also neglect the conformal 
factor. 

• As we already commented, one might be a little worried that the "speed" c depends on the arbitrary but fixed 
scale fjL. But the 3-metric hij also depends on pi and 

c2 h'l = tt" po [m-^]'^ (50) 

is /X- independent. Similarly, {vqY is independent of /i. If you ask physical questions like (for instance) "what 
are the null curves of the metric g^^i/?" they are determined by the equation 

And this equation is completely independent of the arbitrary fixed scale /i. Equivalently 

f dx' ,\ f dx^ ,\ „ , , 
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We should add that in Einstein gravity the spacetime metric is related to the distribution of matter by the non- 
linear Einstein-Hilbert differential equations. In contrast, in the present context, the acoustic metric is related 
to the background wavefunction in a simple algebraic fashion. There are certainly constraints on the acoustic 
metric, but they arise from the generalized nonlinear Schrodinger equation; not from the Einstein equations of 
general relativity. (We belabour this trivial point because we have seen it lead to considerable confusion.) 

Finally we reiterate the main reason the relativity community is interested in these systems: regions where the 
speed of the condensate flow exceeds the speed of sound very closely mimic the key kinematic features of black 
hole physics. We will not repeat the relevant details as they are more than adequately dealt with elsewhere in 
the hterature. (See, for example, |l|, |, |, |, ||, 0, |, |, |lO[ 0, [l|||.) 



Many features of this acoustic metric survive beyond the current quasi-classical approximation, and we now initiate 
a systematic analysis of how much further the model can be pushed. 



VII. LOW-MOMENTUM APPROXIMATION 



The low-momentum approximation is subtly different from the naive quasi-classical approximation. We shall now 
retain the D2 term, but take the approximation that within this quantum potential term the gradients of the 
background are more important than gradients of the fluctuation. (We justify this with the observation that gradients 
of the fluctuation are doubly small, being suppressed by both a factor of fi and a factor of the linearization parameter 
e). Specifically we take 

1 -3/2 r« / +l/2^l 

[^h{Po )\ Pi 



D2PI 



1 _ 



1 -1/2 . , -1/2 ^ 

Po ^h[Po Pi) 



-1/2 



2 

1 _ 



1/2 



Pi 



^hPO (V/9o)2 



Pi 



Pi 



Pi- 



That is, under this low-momentum approximation we can simply replace the operator D2 by the function 



Pi 



Pi 



(53) 
(54) 

(55) 
(56) 



The net result of this approximation is that wherever the quantity n" appears in the naive quasi-classical analysis it 
should be replaced by 



Ap 



^hPQ (Vpo) 



pI 



pI 



This does not affect the background flow velocity vq, but it does modify the propagation speed so that now 



4m 



A/iPo (Vpo) 



Pi 



Pi 



classical 



1 



Afin" 



AftPo (Vpo) 



Pi 



Pi 



(57) 



(58) 



This is a new higher-order correction to the quasi-classical speed of sound we have previously introduced in equation 
( pl| ) ; in this sense is a generalization of the previously known results regarding the propagation of collective excitations 
in BEC. This effect was not contemplated in the Garay et al analysis. (Justifiably so, since they were only interested 
in the geometrical "optics" approximation within the quasi-classical limit.) If, using the speed of sound, we introduce 
a notion of "acoustic Compton wavelength" 

' (59) 



A. 



then 



pc 



= — (60) 

PPO 

So this modification to the speed of sound is seen to be the first order term in a gradient expansion governed by the 
dimensionless parameter 

llVpoll 



A. 



inpo 



(61) 
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VIII. EIKONAL APPROXIMATION 



In contrast to the low-momentum approximation, the eikonal approximation is a high-momentum approximation 
where the phase fluctuation 6i is itself treated as a slowly-varying amplitude times a rapidly varying phase. This 
phase will be taken to be the same for both pi and 9i fluctuations. In fact, if one discards the unphysical possibility 
that the respective phases differ by a time varying quantity, any time-constant difference can be safely reabsorbed in 
the definition of the (complex) amplitudes. 

Specifically, we shall write 



9i{t,x) — Re{Ae exp(— 
pi(i,x) = Re{Ap exp(— 



(62) 
(63) 



As a consequence of our starting assumptions, gradients of the amplitude, and gradients of the background fields, are 
systematically ignored relative to gradients of (p. [Warning: what we are doing here is not quite a "standard" eikonal 
approximation, in the sense that it is not applied directly on the fluctuations of the flcld ^{t,n) but separately on 
their amplitudes and phases pi and 4>i .] We adopt the notation 



Then the operator D2 can be approximated as 

D2 Pi = 



1 _ 



3/2 r. / +1/2m 

[A/i(Po )] Pi 



-1/2 A / -1/2 \ 



1 



2^0 ^ Pi- 



(64) 

(65) 
(66) 
(67) 



A similar result holds for D2 acting on 9i. That is, under the eikonal approximation we effectively replace the operator 
D2 by the Junction 



D2 



1 



For the matrix /'"^ this effectively results in the replacement 

r 





^"{po) + 


pi 


-1 








4^ po. 










vr"(po) + 




-1 


V 






4^ po. 




P 






7r"(po) + 


P PI 






P 


4p Po. 




Po 




h 


^"(po) + 




P 


P 





P 



4p Po 



P 



(As desired, this has the net effect of making f^" a matrix of numbers, not operators.) 
The physical wave equation (^9|) now becomes a nonlinear dispersion relation 

f° ^2 ^ (^0. ^ ^.0) ^ ^ f .j ^ 0^ 

After substituting the approximate D2 into this dispersion relation and rearranging, we see (remember: 



■ UJ 



2 Wq uoki 



pok'^ 



tt" + 



4ppo 



-P 



0. 



That is: 



4ppo 



(68) 

(69) 
(70) 
(71) 
(72) 

(73) 

(74) 
(75) 
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Alternatively 



IV = Vq ki i: 



^vlh± ^/po7r"(A:,[m-i]y/c,) + ( -(k,[m-^Yo k^) ) . (76) 



In the case of an isotropic mass matrix, [m~^]*-' 1/m. It also makes sense to set /i — > m, and hij —> 5ij, in 
which case c really does represent the physical speed of sound. Then we can write the dispersion relation in a more 
illuminating form 



Lo = vi fc,±^c2p+ (^A . (77) 

Notice how the previous anisotropic dispersion relation differs from this isotropic one: in that case there are three 
different physical sound velocities in three orthogonal principal directions. 
At this stage some observations are in order: 

— (1) It is interesting to recognize that the dispersion relation ( |77| ) is exactly in agreement with that found in 1947 
by Bogoliubov |2^ for the collective excitations of a homogeneous Bose gas in the limit T ^ (almost complete 
condensation). In his derivation Bogoliubov applied a diagonalization procedure for the Hamiltonian describing the 
system of bosons. 

— (2) It is easy to see that ([76|), and its isotropic partner ( |77| ) actually interpolates between two different regimes 
depending on the value of the wavelength A = 27r/| |fc| | with respect to the "acoustic Compton wavelength" Xc = h/ {fic). 
(Remember that c is the speed of sound; this is not a standard particle physics Compton wavelength. Furthermore 
is defined using the inverse 3-metric h^^ .) In particular, if we assume vq — (no background velocity), then for large 
wavelengths A ^ Ac one gets a standard phonon dispersion relation « c||A;||. As stressed by Braaten this can be 
related to the fact that the quantum theory we are working with has a U{1) symmetry which is spontaneously broken. 
At low momenta we are just seeing the dispersion relation of the corresponding Goldstone mode. For wavelengths 
A Ac the quasi-particle energy tends to the kinetic energy of an individual gas particle and in fact w ~ h^k'^ /{2m). 

— (3) The dispersion relation ([77| ) exhibits a contribution due to the background flow Vq ki, plus a quartic dispersion 
at high momenta. The group velocity is 



duj 
dki 



vl ± / — k' (78) 



^,2p+(_l.p) 



Dispersion relations of this type (but in most cases with the sign of the quartic term reversed) have been used by Corley 
and Jacobson in analyzing the issue of trans-Planckian modes in the Hawking radiation from general relativistic black 
holes H, 1^, |20). In their analysis the group velocity reverses its sign for large momenta. (Unruh's analysis of this 
problem used a slightly different toy model in which the dispersion relation saturated at high momentum [^.) In 
our case, however, the group velocity grows without bound allowing high-momentum modes to escape from behind 
the "horizon". (Thus the acoustic horizon is not "absolute" in these models, but is instead frequency dependent, a 
phenomenon that is common once non-trivial dispersion is included.) 

This type of "superluminal" dispersion relation has also been analyzed by Corley and Jacobson They found 
that this escape of modes from behind the horizon often leads to self- amplified inabilities in systems possessing both an 
inner horizon as well as an outer horizon, possibly causing them to disappear in an explosion of phonons. This is also 
in partial agreement with the stability analysis performed by Garay et al using the whole Bogoliubov equations. They 
found unstable solutions with the kind of behaviour just mentioned, but they also find stability regions (depending 
on the value of certain configuration parameters). The existence of this stable configurations might be related with 
the specific boundary conditions imposed in their configurations. 

Indeed, with hindsight the fact that the group velocity goes to infinity for large k was pre-ordained: After all, we 
started from the generalized nonlinear Schrodinger equation, and we know what its characteristic curves are. Like the 
diffusion equation the characteristic curves of the Schrodinger equation (linear or nonlinear) move at infinite speed. 
If we then approximate this generalized nonlinear Schrodinger equation in any manner, for instance by linearization, 
we cannot change the characteristic curves: for any well behaved approximation technique, at high frequency and 
momentum we should recover the characteristic curves of the system we started with. However, what we certainly do 
see in this analysis is a suitably large region of momentum space for which the concept of the effective metric both 
makes sense, and leads to finite propagation speed for medium-frequency oscillations. 
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— (4) There is an amusing feature to the (generalized) Bogoliubov dispersion relation which it may be worth making 
explicit: Consider the dispersion relation 

(BEC condensates correspond to mo = 0, we retain this term here for generality. We have made c = fi = 1.) At low 
momenta {k <C TOq) this dispersion relation has the usual non-relativistic limit 

Lo{k)=ma + - + 0{k^). (80) 

At intermediate momenta (mo k <^ ttLoo) this dispersion relation has an approximately relativistic form. Finally 
at large momenta (fc 3> moo) the dispersion relation again (perhaps surprisingly) recovers a non-relativistic form 

cc'(fc) = -^+moo + 0(fc-2). (81) 

This serves to drive home in a particularly simple way the point that observing a Lorentz invariant spectrum does 
not guarantee that the underlying physics is Lorentz invariant. Indeed the entire programme of searching for analog 
models of general relativity generically seeks to take some simple (physically accessible and typically non-relativistic) 
model and determine if it nevertheless hides within it some useful approximation to Lorentzian geometry. 



IX. SUMMARY AND DISCUSSION 



In seeking to see how far we might be able to push the BEC system as an analog model for general relativity we 
have encountered a number of intriguing features: First, we have shown that the existence of a regime in which phase 
perturbations of the wave function of the BEC (or quasi-BEC) behave as though coupled to an "effective Lorentzian 
metric" is a generic feature, independent of the explicit form of the non-linear terms in the Schrodinger equation. The 
only exception comes about when the forces exerted between atoms are attractive. In this case the equation of motion 
for the phase perturbations are no longer hyperbolic, so the whole notion of wave disappears for these systems. (In 
GR language, this corresponds to a Euclidean-signature metric.) 

Second, we have seen that in contrast to the isotropic acoustic systems considered to date, mimicking a generic 
gravitational field is not a priori implausible though it would technically be very challenging, relying as it does on 
the direct introduction of anisotropics into the generalized nonlinear Schrodinger equation via a position-dependent 
3-tensor effective mass. 

We have also explicitly seen how the whole notion of "effective metric" in these BEC systems is intrinsically an 
approximation valid for certain ranges of frequency and wavenumber — this should not really be a surprise since even 
for normal acoustic systems eventually the atomic nature of matter provides a natural cutoff. In the BEC system 
we have seen that the acoustic Compton wavelength plays a similar role: wavelengths long compared to the acoustic 
Compton wavelength see a Lorentzian "effective metric" while wavelengths short compared to the acoustic Compton 
wavelength probe the "high energy" physics (which in this situation is the non-relativistic Schrodinger equation). 

Finally, we mention that from the sound velocities typically encountered in BEC systems we can make a crude esti- 
mate of the Hawking temperature to be expected in these systems. Using the standard estimate based on dimensional 
analysis 

and choosing a value of i? w 10 /i for the size of the acoustic black hole, we would have T « 10^^ K. While extremely 
small this temperature is only three orders of magnitude less than that of the BEC itself. Furthermore, as argued 
in | p2| this order of magnitude estimate for the Hawking temperature is often misleadingly low — and this fact, in 
addition to the intrinsically interesting features of the effective geometry approach, makes further investigation of 
these BEC systems worthwhile — both for the condensed matter and relativity communities. 
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